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We investigate the effect of spacetime backreaction on tfie upper critical magnetic field for s-wave 
holographic superconductors by using the matching method. The backreaction of the constant external 
magnetic field and the electric field to the background geometry leads to a dyonic black hole solution. 
The magnetic fields as well as the electric fields acting as gravitational sources tend to depress the 
critical temperature of the superconductor. We derive the analytical expression for the upper critical 
■ magnetic field up to 0{k^) order and find that backreaction makes the upper critical magnetic field 

stronger. The result is consistent with the previous numerical and analytical results. 
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\\. Introduction 



The gauge/gravity duality Ull-ISll] as the most fruitful idea stemming from string theory, 
has been proved to be a powerful tool for studying the strongly coupled systems in field 
theory. By using a dual classical gravity description, we can effectively calculate correla- 
tion functions in a strongly interacting field theory. Recently, a superconducting phase was 
established with the help of black hole physics in higher dimensional spacetime ©-[T])]. 
^ , Counting on the numerical calculations, the critical temperature was calculated with 

and without the backreaction for various conditions Il8ll-l24l)1. The behavior of holographic 
^ ■ superconductors in the presence of an external magnetic field has been widely studied in 

^ , the probe limit Il25] )-i35])1. The analytical calculation is useful for gaining insight into the 

^ \ strong interacting system. If the problem can be solved analytically, however vaguely, one 

^ ■ can usually gain some insight. As an analytical approach for deriving the upper critical 

! magnetic field, an expression was found in the probe limit by extending the matching method 

Xf^ \ first proposed in 19)] to the magnetic case Il32l )1. which is shown to be consistent with the 

O ■ Ginzburg-Landau theory. 

Most of previous studies on the holographic superconductors focus on the probe limit 
neglecting the backreaction of matter field on the spacetime. The probe limit corresponds 
to the case the electric charge g — oo or the Newton constant approaches zero. The back- 
^ \ reaction of the spacetime becomes important in the case away from the probe limit. At 

a lower temperature, the black hole becomes hairy and the phase diagram might be mod- 
ified. Recently, an analytical calculation on the critical temperature of the Gauss-Bonnet 
holographic superconductors with backreaction was presented in 121 1) 1 and confirmed the 
numerical results that backreaction makes condensation harder |[T3T l.[T4]|.[T9]l.l20l)1 . 

Considering the above facts, it would be of great interest to explore the behavior of 
the upper critical magnetic field for holographic superconductors in the presence of the 
backreaction. In this work, we will first consider the effect of the spacetime backreaction to 
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s-wave holographic superconductors without the magnetic field. Different from the probe 
limit case, the backreaction of spacetime actually leads to a charged black hole solution 
in AdS space at the leading order. We will compute the critical temperature analytically 
by using this charged black hole metric through the matching method. Secondly we will 
study the properties of holographic superconductors in the presence of external magnetic 
field. When we turn on the external magnetic field, the resulting background geometry 
becomes the dyonic black hole solution in AdS space to the zeroth order. The analytical 
investigation on the effect of the spacetime backreaction to the upper critical magnetic field 
has not been carried out so far. Therefore, the contents in this paper will be greatly different 
from the probe limit case as we consider the spacetime backreaction. Note that in both 
cases, the small backreaction approximation shall be used to obtain an analytical result. In 
order to compare the analytic study with numerical results, we will also carry on numerical 
computation. 

The organization of the paper is as follows: We first consider the effects of backreaction 
on 2 + 1-dimensional s-wave holographic superconductors in section 2. Without the mag- 
netic field, the critical temperature with backreaction will be derived first. Then we continue 
the calculation to the strong external magnetic field case and find an analytical expression 
for the backreaction on the upper critical magnetic field. From the Einstein equation, we 
know that the presence of charge and magnetism in 4-dimensional spacetime yields a dy- 
onic black hole solution. The critical temperature may influenced by the backreaction of the 
magnetism. We will compai^e the analytic and numerical results. We present the conclusion 
will be presented in the last section. 

§2. (2 + 1) -dimensional s-wave holographic superconductors 

In this section, we first investigate the backreaction of electric field on superconduc- 
tivity and derive the phase transition temperature Tc in this case. After that, we turn to the 
backreaction of the external magnetic field and calculate the critical magnetic field. 

2.1. Critical temperature with backreaction in the absence of magnetic fields 

We begin with a charged, complex scalar field into the 4-dimensional Einstein-Maxwell 
action with a negative cosmological constant 

S = j d\^gi^R -2A- - \d^^l^ - iqA^i^l^ - mVl'}, (2-1) 

where G4 is the 4-dimensional Newton constant, the cosmological constant A = and 
F/^iu = dfj,Ai, — dyA^. The hairy black hole solution is assumed to take the following metric 
ansatz 

ds" = -f{r)e-^^^Ut' + 7^ + j^idx^ + dy^), (2.2) 

together with 



^^ = ((/<(r), 0,0,0), V = V'(r-). 



(2.3) 
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The Hawking temperature, which will be interpreted as the temperature of the holographic 
superconductors, is given by 



T 



-x(r)/2 



(2.4) 



r=r+ 



where a prime denotes derivative with respect to r and r+ is the black hole horizon defined 
by f{r^) = 0. ipir) can be taken to be real by using the U{1) transformation. The gauge 
and scalar equations become 



2j,2 



f y' 2 
< + (y-| + -)V'' + 



0, 



P 



m . 



The tt and rr components of the background Einstein equations yield 



/ 



3r 2 



L'2 



+ mV + /( + 



1.12 



0. 



0, 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



When the Hawking temperature is above a critical temperature, T > Tc the solution is the 
well-known AdS-Reissner-Nordstrom (RNAdS) black holes 



/ 



r 



P 



+ 



2r+ 



+ 



X = ip = 0, (/) = p 



1 1 



(2.9) 



where = 87rG4. At the critical temperature T = the coupling of the scalar- to gauge 
field induces an effective negative mass term for the scalar field,the RNAdS solution thus 
becomes unstable against perturbation of the scalar field. At the asymptotic AdS boundary 
(r — )• oo), the scalar and the Maxwell fields behave as 



> 



4> = 



+ ... 



(2.10) 



where fi and p are interpreted as the chemical potential and charge density of the dual field 
theory on the boundary. According to the gauge/gravity duality, < Oa± > represents the 
expectation value of the operator C'z\± dual to the charged scalar field V'- The exponent A± 
is determined by the mass as Z\± = I ± iV9 + 4m2. Note that for ip both of the falloffs 
are normalizable and we choose the boundary condition that either < C'/\_ > or < Oa+ > 
is vanishing. We will impose that p is fixed and take < >= as in O]. Moreover, we 
will consider the values of which must satisfy the Breitenlohner-Freedman (BF) bound 



^t-'bf — < '^BF + 1 ^^^^ ''^'bf ~ ~{d — l)/4 ll37] )1 for the dimensionality of the 



^BF 



spacetime d = 4 in the following analysis. 
After introducing the new coordinate z 



-/' + 



/ 3r2 



Pz3 



, the equations of motion become 

4 g2^2^2gx 



0, (2.11) 
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2 r ^4 



.2^ ' p 



0, (2.12) 



^" + ^X^'- J^^ = O, (2.13) 



<-(x'-f)^'4(^-f )^^0, (2.14) 



where the prime ' denotes a derivative with respect to z. One may find that the transformation 
(p ^ (p/q and ip ^ ip/q does not change the form of the Maxwell and the scalar equations, 
but the gravitational coupling of the Einstein equation changes — >• jq^. The probe 
limit studied in 17])] corresponds to the limit (7 — )• oo in which the matter sources drop out of 
the Einstein equations. The hairy black hole solution requires to go beyond the probe Umit. 
In l?])], it was suggested to take finite q by setting = 1. Recently, the author in ||2T1 )1 
proposed to keep 2k;^ finite with setting q = \ instead. We will take the latter choice. 

In the neighborhood of the critical temperature Tc,we can choose the order parameter 
as an expansion parameter because it is small valued 

e=<OA^> ■ (2.15) 

We find that given the structure of our equations of motion, only the even orders of e in the 
gauge field and gravitational field, and odd orders of e in the scalar field appear here. That 
is to say, we can expand the scalar field V', the gauge field as a series in e as 

(/> = (/>0 + 6^02+ 6^4 + •■• (2.16) 

V' = e'0i + e^V3 + e^V'5 + --- (2-17) 

Let us expand the background geometry line elements /(z) and x(z) around the AdS- 
Reissner-Nordstrom solution 

/ = /o + e'/2 + eV4 + ... (2.18) 

X = e'x2 + e'x4 + - (2.19) 
The chemical potential \i should also expanded as 

/i = /io + £^(5/^2, (2.20) 

where 5/^2 is positive. Therefore, near the phase transition, the order parameter as a function 
of the chemical potential, has the form 

It is clear that that when \i approaches /^o, the order parameter e approaches zero. The phase 
transition occurs at the critical value [ic = Ho- Note that the critical exponent 1/2 is the 
universal result from the Ginzburg-Landau mean field theory. The equation of motion for 
is solved at zeroth order hy (po = /Uo(l — -2) and this gives a relation p = /xor+. So, to zeroth 
order the equation for / is solved as 



Uz) = ^(l-z)(l + z + z^- ^z^ ] . (2.22) 
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Now the horizon locates at z = 1. We will see that the critical temperature with spacetime 
backi^eaction can be determined by solving the equation of motion for i/j to the first order. 

At first order, we need solve the equation for ipi by the matching method. The boundary 
condition and regularity at the horizon requires 



V^'(i) 



2 2 



/o(l) 



V'i(i). 



In the asymptotic AdS region, ipi behaves like 

Now let us expand tpi in a Taylor series near the horizon 

V^i = Ml) - + ^V/i'(i)(i - zf + ... 

From (12. 141 i. we obtain the second derivative of tpi at the horizon 



<(1) 



; 4 + 



2 2 

m r± 



m 

/o(l) /o(l) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



Using (12.231 ) and (12.261) . we find the approximate solution near the horizon 



V'i(z) = ^i(l)- 



2 2 



2 2 



4/^(1) V m m)) 4 aiy 



\21 



Mm-zf+. 



{121) 

In order to determine and C+, we match the solution (12.24b and (12.251 ) smoothly at 

Zm- We find that 



z^+C+ = V'i(l) 



2 2 



2 2 



4/^(1) 



4 + 



/o(l) 



2 2 



A+z^+-'C^ 



4 

2 2 

r±m , , 

r2 



(2.28) 



2 2 / 



4/, 



^i(l)(l -z^). 



4 + 



/o(l) /o(l) 



4 m)\ 

Solving the above equation, we obtain the expression for C+ in terms of "01(1) 



1 



1 - Zm rlm'^ \ 



/( 



Ml)- 



(2.29) 



(2.30) 



2zm + {l-z„^)A+ 

Substituting the above equation back into ( 12.291 ). we get a non-trivial relation provided 

V'l(l) 7^0, 



2A. 



2Zm + (1 - Zm)A^ 



(1 - Zra)A^ 



2Zm + (1 — Zjn)A^ 



+ (3 - 2z„ 



n 2 2 



m 



(1 - Zm)rlm^ /^(l) 1 - (1 - z^)rl <p'^{lf 



m)' 



+ 



my 



my 



0. (2.31) 
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Note that /^(l) 



2ri 



these relations back into (12.311 ). we obtain an equation for fiQ 



6+^^ ) and</.^(l) = -^0- Plugging 



4;4 



kH 



A, 



2r% 2Zrr, + (1 - Zm)A+ 



4 



2r2 



1 + 2k' 



- z„ 



2H(2z„ + (l-z^)Z\+) 



2;4 



r4(2z^ + (l-z^)Z\+) 
4;4 



^0 



+37712^2 (2 - Zm) + 



ml 



18 + 3m2p(l-z„) 



0. 



(2.32) 



z„(Z\+ - 2) - Z\ 

The main idea of ||2T1 )1 is to work in the small backi-eaction approximation <^ 1 together 
with the matching method so that all the functions can be expanded by and the term 
in the above equation can be neglected. In this sense, fiQ is solved as 



fJ'O 



■2m. 
2/2/ 



4)4 



ml 



18 + 3m2/2(i 



+ 



Zm{A+ - 2) 

3A~ 



A. 



1/2 



2 

2^2 



2Zm + (1 - ^^m)^- 



P{1 - Z, 



mH\\ - Zm)A^ 

A{2Zm + (1 - Zm)A^ 



(2.33) 



Without the tenii, the expression for /io can be reduced to the result of the probe limit 
case. The term in the above equation is positive, which means that //q increase. By 
further using the relation /^o = — , we find an expression for r+: 



r+= p 



a/4 



m 



2,mH\2-Zm) + ^—{l-z, 



+ 



18 + 3m2/2(l - 
■ Zm{A+-2) - A+ 
3A+ 



■A, 



1/4 

'iZffi 



1 + 



2k" 



l\l 



m^l 



2;2 



8(2z„, + (1 - Zn,)A^ 



2{2Zm + (1 - 

The Hawking temperature is given by 

T 



2r2 
m I 



2rl 



(2.34) 



(2.35) 



When /io = /ic, the above Hawking temperature reaches the critical point Tc where the order 
parameter approaches zero. From (I2.33t and (I2.35t . we obtain the critical temperature 



2,m^l\2-Zm)- 



m I 



1 1 Znri 



l^ + 3mH\l - z. 



Zm{A^ 



A, 



Together with (12.341 ). we write the critical temperature in a form as 

2^2 



Tc = Ti{l 



,2-^2) 



-A, 
(2.36) 

(2.37) 
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where 



^ 3pV2/i 



47r/ V 2 



37n^f(Zra-2) + ^-—(l 



ml 



4;4 



18 + 3m2/2(l_2^)^ - 



,(Z\+-2)-Z\. 



4 



(2.38) 



36Z\+ + m2/2 ^ m^PA+ 



8(1 - Z„^)[2Z^ + (1 - Z^)A+] 2[2Zm + (1 - 

(12 - 7zmWl^ rriH^ 

+ (2.39) 



8(1 - z^) 12 

It is easy to check that when m'^P = —2, = 1/2 and thus Z\+ = 2, we have Ti = 
— — , the exact resuh obtained in ©l for (2 + 1) -dimensional superconductors and 

T2 = |. This result is also in good agreement with the numerical result by choosing a 
proper matching point Zm 17])]. We also find that the corrections due to the backi^eaction T2 
is positive for arbitrary Zm in the region (0 < < 1)- Therefore, we confirm the result 
found in |[T9l )-l2TI)1 that the backreaction makes condensation harder. The reason for the 
decreasing of the critical temperature can be understood from the relation Tc oc 1/ i_Iq 
The value of fiQ increases due to the gravitational backreaction and thus Tc decreases. 

2.2. The upper critical magnetic field with backreaction 

In this section, we will explore the effects of the backreaction on the external critical 
magnetic field. In the neighborhood of the upper critical magnetic field Bc2, the scalar field 
ip is small and can be regarded as a perturbation. The scalar field becomes a function 
of the bulk coordinate z and the boundary coordinates (x, y) simultaneously because of the 
presence of the magnetic field. According to the AdS/CFT correspondence, if the scalar field 
-0 ~ X{x, y)R{z), the vacuum expectation values < O >oc X{x, y)R{z) at the asymptotic 
AdS boundary (i.e. z — )• 0) |[25l) .[30l)1. We can simply write < O >oc R{z) by dropping the 
overall factor X{x, y). So, to the leading order, it is consistent to set the ansatz 

At = Mz), = 0, Ay = Bc2X. (2.40) 

Note that the applied external magnetic is constant and homogenous. Considering the fact 
that an external magnetic field is included, one may wonder whether such a constant external 
magnetic field could backreact on the bulk gravity or not. We may need to consider the 
effects of the spatial component of the gauge field in the superconducting phase and assume 
the gauge field behaves as 

A = 4)Q{z)dt + h{z)dx. (2.41) 

In other words, we need solve the bulk gravity equation for h{z) and the resulted metric 
is anisotropic, following this line, we may obtain a kind of dyonic black hole solution, 
which includes chai^ge and magnetism. However, we notice that several authors have already 
discussed such conditions in |[36l) 1. In these works, h{z) is interpreted as the vector hair of 
the black hole. At the AdS boundary, h{z) behaves as 

h{z) = a - iz + ... (2.42) 
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According the AdS/CFT correspondence, ^ is the dual current density and a is the dual 
cunent source of the holographic supeifluid. Of course, it is not proper to regai^d ^ as a 
homogenous applied magnetic field. Actually, when we discuss the vortex structure of the 
holographic superconductors, in general we should consider 



ipi = ipi{x,y,z), At = (p{x,y,z), = A^{x,y, z), Ay 



1 -^c/ 



Ay{x,y,z). (2.43) 
A^{g, z) as well as 



or simply in the polar coordinate V'l = i^i{g, z), At = 4>o{q, z), 
the boundary condition that At{z = 1) = {) and Aip{z = 1) regular. 

In the presence of external magnetic field, not only the matter fields but also the space- 
time metric should depend on the coordinates {z,x,y). The background static metric may 
have the fonii 



goo{z,x,y)dt^ + g;,Az, x,y)dz^ + g^^{z,x,y)dx^ + gyy{z, x,y)dy^ 



In this case, we need solve the Einstein equations 



R 



ilM) 



(2.45) 



where 



together with the Klein-Gordon equation 

1 



and the Maxwell equation 



-g 



-g 



--dx 



gg^'D.ij 



-gg^^gf'F, 



(2.46) 



(2.47) 



(2.48) 



where we have defined D^ip = d^tp — iA^ip and Jo- = i[ip*D^ip — 'tp{Da-tp)*]. In this 
case, we have three coupled nonlinear partial differential equations involving the metric 
components, scalar field ip, the scalar potential At and vector potential A in which analytic 
study becomes very difficult to do. Note that we can expand the background geometry in 
series of e 



.(0) 



2„(2) , ^4^(4) 



g^.u = g'fTJ + e-grj + e^g';^ + ■■ 



(2.49) 



To solve these equations analytically we will follow the logic as shown in table Jl In the ab- 
sence of the external magnetic field, the backreaction of the electric field to the background 
geometry leads to the RNAdS black hole solution at the zeroth order. At the linear order, 
the metric receives no corrections from matter fields and we need only solve the equation 
of motion for Tpi at this moment. As we have done from (12.141 ) to (I2.371 i. all the equations 
depends only on the radial coordinate z. We obtained the critical temperature with back- 
reaction. When we turn on the external magnetic field, the background spacetime changes 
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because of the presence of Bc2- We can still expand ip, At and in series of e. At the 
leading order, the matter field (po and A^^ result in a dyonic black hole solution in AdS 
space. By solving ■ipi(xi,z) (xi = x, y) at next to leading order, we should obtain the ex- 
pression for the upper critical magnetic field. The above arguments are actually the logic of 
the calculation of the whole paper. 





Vanishing magnetic field 


External magnetic field 




e'i}i{z) + e-^'ili3 + ... 


6^1(2:, Xi) + e^ili^iz^Xi) + ... 


At = 


e>oW + e"<^2(2) + ... 


e"<^o(z) + e''02(z,2;i) + ... 









gt^<^ = 


e^QRNAdS + <^^g\tu + ■■■ 


t y dyonic \ fc y/^/v ■ • ■ 



Table I. Logic of the analytic calculation. 

After justify the usage of (I2.401 i. we can then solve the equations of motion order by 
order. The black hole carries both electric and magnetic charge and the bulk Maxwell field 
yields 

A = Bc2xdy + (pdt (2.50) 

At the zeroth order 0{e^),vJt solve the Einstein equation and the line elements of the dyonic 
black hole metric ai^e given by ||38] )1 



fodt' + -fp^{dx' + dy') + ^dz\ 

_ _p_ 



(2.51) 
(2.52) 



where /o = ^|2- ^1 — ^1 + z + — ^^,2° z^ — ^^^^^z'^^ . The Hawking temperature 



at the event horizon is evaluated as 



T 



2ri 



c2 



2ri 



(2.53) 



To the linear order, the equation of motion for ipi has its new form 



m I ip 



dl + {dy - iB,2xf 



(2.54) 



where the prime denotes derivative with respect to z. We use sepai^ation of variables 

^1 = e"'yyXn{x)Rn{z), (2.55) 
and obtain the equation of a two dimensional harmonic oscillator and a equation for R{z) 



+ (ky - B,2xfXn{x) = \nB,2Xn{x), 



foR'n + fo^'n + 



fo 



m ]Rn 



XnBp.oP 



(2.56) 
(2.57) 
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where A„ = 2n + 1 is the eigenvalue of the harmonic oscillator equation, n = 0, 1, 2, ... 
denotes the Landau energy level and the prime in (12.561 1 and (12.571 ) denote derivative with 
respect to x and z, respectively. The equation (I2.56l )is solved by the Hermite polynomials 

(Bc2^-fcy)^ 

Xn{x) = e 2«c2 Hn{x). (2.58) 

Let us choose the lowest mode n = in what follows, which is the first to condensate and 
is the most stable solution after condensation. Actually, the Arikosov vortex lattice is given 
by a superposition of the lowest energy solutions 

V^i = i?o(z)5^c,e*'=^^Xo(x), (2.59) 
j 

where cj are coefficients that determine the structure of the vortex lattice. 

Now we are going to solve (12.571 ) by exploring the matching method and find the cor- 
rection to the upper critical magnetic field away from the probe limit. Again regularity at 
the horizon requires 

^o(l) = ^^o(^) + ^^-^^^ 
The behavior of Rq at the asymptotic AdS boundary is given by 

Ro{z) = C+z^+. (2.61) 

The scalar potential (f)Q satisfies the boundary condition at the asymptotic AdS region (j)o{z) = 
— :f^z and vanishes at the horizon (/iq = 0, as z — )• 1. In the strong field limit, the scalar 

field ^|) is almost vanishing and we can drop out the jV'P term in the right hand side of equa- 
tion (12.131 ). One may find that (j)o{z) = ^(1 — z) is a solution that satisfies (12.131 ) and the 
corresponding boundary conditions Il30] l1. 

In the presence of the external magnetic field, the Taylor expansion of i?o near the 
horizon still goes as 

Ro{z) = i?o(l) - R'oim -z) + - zf + ... (2.62) 

From (I2.57I I. we know that near z = 1, R"{1) is expressed as 

(2.63) 

Putting the expressions for Rq{1) and -Ro(l) into (12.621) . we obtain 



Ro{z) = i?o(l) - ( ^ + ^ )i?o(l)(l -z) + 



/^(l) J 2/^(1) 4 /^(1)2 



4/^(1) V m 

Ro{l){l - zf + ... (2.64) 
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We connect the two solutions (12.611 ) and (12.641 ) at a intemiediate point Zm smoothly and thus 
find that 



1 1 



+ 



Roil) ■ 



+ 



Ro{l){l - z„,) + 
rim? 



2/^(1) 4 /^(l)^ 



4/^(1) 

Ro{l){l-Zmf, 



(2.65) 



^W + B,2P 



+ 



B,2l^ r2</>;(l)^ 



lm^ + B,2l^ ( /^(l) rlw? + B,2l^ 



4/^(1) 



4 + 



/o(l) 



/o(l) 



2/^(1) 4 
From the above equations, we find that 

22m — Zi-u 



C4 



2Zm, + (1 - Zra)A. 



i?o(l)(l-Zm). 

_ 1 - Zm + Bc2l'^ 

2 /^(l) 



(2.66) 



(2.67) 



Substituting the above relation back into (12.661 ). we get a non-trivial expression 



2Zm + (1 - 2;m)Z\H 



(1 - 2m)^H 



22;^ + (1 - ^;m)^H 



+ (3 - 2Zm) 



/o(l) 



2\2 



(1 - Zm)Km^ + B,2n /^'(l) 1 - Zm + B,2n 



+ 



+ ^77Tt(1 - Z„,) 



m' ■ 2 

(1 - zm)ri 0'o(l)' 



0. 



(2.68) 



/^(l)^ 2 /^(1)2 

When we turn off the magnetic field Bc2 = 0, (12.681 ) returns to (I2.311 l. In the presence of 
the magnetic field, both the charge and the magnetic field can backi^eact on the black hole. 
The critical temperature should receive further connections from the magnetic field. The 
difference between (12.681 ) and (12.321 ) comes from the Bc2 related terms, which goes as 

(Zm - 



2Z\. 



2Zm (1 ^m) 

1 - Zr, 



■'c2 



Zm{A+ - 2) - Z\h 



2 



2rl 



Wc2n + 
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^0- (2-69) 



We obtain a relation between Bc2 and r+ by using (12.331) and 
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The critical temperature dropped because of the magnetic field 
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), 



1/2, q = 1, 
(2.70) 

(2.71) 
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where Ti 



3v^ 



This reflects the fact that condensation becomes even harder when 
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one turns on the external magnetic field. 

Note that (12.701 ) is not enough to determine the relation among the upper critical mag- 
netic field, the system temperature T and the critical temperature Tc. Considering the values 
of /o(l), /o (1) and 4>q{1) and solving (12.68b to the first order of k^, we get 



n 



1^0 



K 



1 



R2 /6 



1 - Zm nr[ 

+2Z\+ + m2/2(2(8 + ^Zm{A+ - 2) - 8Z\+) + ^A+)/rl^ [zm{A+ - 2) - Z\+]~^^ | 
with 

6m2/2 ( z^(4 + Zra{A+ - 2) - 4Z\+) + 3Z\ 



12Z\+ + 2B,2l* zmi'i + z^{A+ - 2) - 3Z\+) 



(2.72) 



+ 



+ m*r(z™ - 1) + 



2) - A, 

12B,2(^m + ^+(l - .Zm)) 

r2(z^(Z\+-2)-Zl+) 



z„(Z\+ - 2) - Z\+ 



(2.73) 



Combining the above equation with /ip = we can obtain a relation between Bc2 and r 



We then substitute (12.331 ) and (12.701 ) into the Hawking temperature T = 



2rl 



and now the Hawking temperature plays the role of the critical temperature in 



the presence of magnetic fields. In order to have a clear picture, by choosing m?l'^ = — 2, 
Z\_l_ = 2 and z^. = 1/2 and further using the relation between r+ and T from the new 
Hawking temperature, we find that the upper critical magnetic field Bc2 yields 
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Bi + B2K^ 
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(2.74) 



In this case, the charge density can be evaluated from (I2.71l i. that is p = ^^^^^ ^1 



2«2 



T2 ) . The upper critical magnetic field Bc2 in series of k. can be expressed as 
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Fig. 1. (color online) The coefficient of the term of the upper critical magnetic field as a function of the 
temperature T/Tc. We set Tc = 1 here 
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Fig. 2. (color online) For fixed external magnetic field, the phase transition temperature decreases if be- 
comes larger. We set Tc — 1 here. 

It is worth noting that Tc means the critical temperature without magnetic fields and 
gravitational backreaction. The result (12.751 1 also implies that it is only applicable near 
the critical temperature Tc because the term will be divergent in the low temperature 
Mmit. One may find that when = 0, the result exactly agrees with ||32] |1. which is also 

consistent with the Ginzburg-Landau theory where Bc2 oc ^1 — T/Tc^ . We also find that 

the coefficient of the term is positive for the system temperature T (see Fig. [T]). This 
result indicates that the effects of the backreaction enhance the value of the upper critical 
magnetic field. The increasing of the magnetic field Bc2 can be explained from the relation 
that Bc2 oc fiQ for fixed value of r_|_ |[32l) 1. Therefore, if the value of fiQ becomes larger, then 
Bc2 increases. However, this does not mean condensation become easy in the presence of 
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the magnetic field. We can see from Fig. |2]that for fixed magnetic field, the phase transition 
temperature goes down as increases. 

2.3. Numerical results 
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Table II. The critical temperature T drops as n increases in the absence of the magnetic field. 
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Fig. 3. (color online) Left: The external magnetic field as a function of T/Tc at = (Green) and k'^ = 0.01 
(Blue). Right: The phase transition temperature decreases if increases in the case that Bc2 = 2. In both 
cases, we choose Tc = 1. 



For completeness of our study, we carry on numerical computation in this subsection. 
We first solve the equations of motion (12.111 ) to (l2.14l ) in the absence of the external magnetic 
field and from which we can obtain the critical temperature and the phase diagram. The 
properties of holographic superconductors without magnetic fields away from the probe 
Umit were studied numerically in llT3l l1 by setting = 1 and finite q. We work in the case 
q = I but finite instead and set r+ = 1 and Z = 1 in the numerical computation. The 
critical temperature T as a function of the backreaction is shown in Table |lll It is clear 
that the critical temperature T drops as increases, which is in consistent with llT3] ).[T9l)1. 

We then consider the behavior of the external magnetic field numerically to the linear 
level by solving equation (12.54b . In FigjS] (left), we find that the external magnetic field 
drops in different ways for = case and = 0.01 case. This is in consistent with the 
analytic calculation in the range T ~ Tc. That is to say, although the critical temperature is 
significantly suppressed by a non-zero k^, the upper bound of Bc2 become larger. In Figj3] 
(right), we also plot the phase diagram of the critical temperature against the gravitational 
backreaction. When we fix the magnetic field, the phase transition temperature is depressed 
as increases, which is also comparable with the analytic results at qualitative level since 
the analytic method closely depends on the matching point. Note that the numerical results 
presented here can be regarded as a side note because we mainly deal with analytical calcu- 
lation in this paper. A more general and thorough numerical computation in the presence of 
external magnetic field with backreaction is called for in the future. 
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§3. Conclusion 

In this paper, we have investigated the effect of backreaction to the upper critical mag- 
netic field of (2 + 1) -dimensional holographic superconductors in Einstein gravity by using 
the analytical method developed in |[9l).l2T])1. As a consistent check, we have derived the 
critical temperature with backreaction in four-dimensional Einstein gravity and confiiTned 
the numerical result given in |[T3] )1 that backreaction makes the condensation harder to form. 
We have obtained the spatially dependent condensate solutions in the presence of the mag- 
netism. The coefficient of spacetime backreaction on the upper critical magnetic field is 
positive for Einstein gravity, which indicates that the magnetic field becomes strong with 
respect to the backreaction in consistent with reference Il39b l . We have also shown the 
corresponding numerical results for each case. 

We can see that the spacetime backreaction presents us an interesting property of holo- 
graphic superconductors: While the backreaction causes the depression of the critical tem- 
perature, it can enhance the upper critical magnetic field. The upper critical field Bc2 is an 
important parameter because it determines the value of the coherence length and strongly 
affects the critical current density J^. The improvement in Bc2 has been the main research 
topic for some experiments. In this paper, we work in the small backreaction limit (i.e. 

<C 1). So if we regard the backreaction as a factor of "doping" in holographic super- 
conductors, then we may find that the backreaction plays the same role as carbon doping 
in MgB2 reported in recent experiments [40 )1 : It results in the depression in T^, while the 
Bc2 performance is improved. Otherwise, we can treat the probe limit approximation as the 
"effective doping": comparing with the superconducting properties with backreaction, the 
probe limit approximation improves the critical temperature but reduces the upper critical 
magnetic field. In microscopic models of high temperature superconductors, the interac- 
tion between doping and electrons contributes a potential term in the Hamiltonian and the 
self energy of the superconducting quasi-particles will be changed. The self-energy can be 
calculated by using the Green function and the collection to the critical temperature can be 
read off from the Green function Bil l. For holographic superconductors, the effective mass 
term is changed with the variation of k^. The extension of this work to the five-dimensional 
Gauss-Bonnet gravity case would be interesting Il42l) -l47ll1. But since the resulting metric is 
anisotropic and analytic calculation become very difficult and involving, we would like to 
leave it to the future publication by using numerical calculation. 
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